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ABSTRACT 



Context. This paper reports on an XMM-Newton observation of the Seyfert 1 Galaxy Mrk 335 performed as a target of opportunity 
when the source was in an unusually low flux state. 

Aims. The low level of continuum emission unveiled an underlying line-rich soft X-ray spectrum, which can be studied with the 
Reflection Grating Spectrometer. 

Methods. The emission features were analysed at high resolution with unprecedented detail for this source. Line ratio diagnostics from 
H-like and He-like ions indicate that the line emission arises in X-ray photoionised plasma. Extensive simulations were performed 
with the CLOUDY photoionisation code. The physical properties of the line emitting material were derived from the comparison of 
the expected and observed line intensities. 

Results. Because of the degeneracy in the ionisation parameter, a number of different solutions for the electron density and column 
density of the gas are consistent with the spectral diagnostics. This prevents us from uniquely determining the properties of the plasma; 
however, the location(s) of the X-ray photoionised gas can be constrained to within 0.06 pc (outer boundary). This limit places the 
X-ray line emitting gas at the inner edge of the material where the broad optical lines are produced (broad line region). 

Key words, galaxies: active - galaxies: Seyfert - galaxies: individual: Mrk 335 - X-rays: galaxies 



1. Introduction 

The soft X-ray spectrum of Seyfert galaxies has represented for 
astronomers a region of great interest and, at the same time, a 
sort of challenge. With the advent of high-resolution spectrome- 
ters onboard XMM-Newton and Chandra, it has been possible 
to reveal the nature of emission in this spectral band at least 
in obscured sources (Guainazzi & Bianchi 2007 and references 
therein). Many authors concur to indicate that the origin of the 
soft X-ray emission in these sources is caused by photoionisa- 
tion of extended circumnuclear gas on kpc scale irradiated by 
the active nucleus which produces He and H-like transitions of 
heavy elements and L-shell transitions from Fe (e.g. NGC 1068, 
Kinkhabwalaet al. 2002; Mrk 3, Bianchi et al. 2005, NGC 4151, 
Armentrout et al. 2007, to cite the most outstanding cases). 

On the other hand, the scenario for unobscured AGN is not 
so definite. The lack of line-of-sight obscuring material in this 
type of sources, allows the whole nuclear radiation to escape 
unblocked to the observer, so that the soft X-ray spectrum is 
dominated by a smooth excess of continuum rather than features 
from reprocessed material (Piconcelli et al. 2005, Crummy et al. 
2006). 

In about half of the Seyfert 1 galaxies (Crenshaw et al. 2003, 
Blustin et al. 2005), the presence of intervening ionised gas 
is directly observable as a series of absorption features, which 
present velocity shifts when originating in outflowing or inflow- 
ing gas. For the sources studied with the longest integration 
times, it has been possible to place some constraints on the elec- 
tron density of the gas and, consequently, on the distance of 
these warm absorbers, mainly through variability studies. The 
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emerging picture is far from being a homogeneous description 
of the circumnuclear ionised gas. Among many cases, we recall 
that the distance of the absorber(s) in NGC 3783 was estimated 
to be 1-3 pc from the nucleus (Netzer et al. 2003 and Behar et 
al. 2003), whereas it was found to be consistent with subparsec 
scale in NGC 3516 (Netzer et al. 2002) and it was estimated to 
be within 4 light days from the nucleus in NGC 405 1 (Krongold 
et al. 2007). 

Mrk 335 is a Seyfert 1 galaxy (z=0.026) with a long X- 
ray history, since it has been observed by almost all X-ray ob- 
servatories in the past. EXOSAT (Turner & Pounds, 1989) and 
BBXRT (Turner et al. 1993a) revealed the presence of an excess 
of emission in the soft X-rays band, later confirmed in ASCA 
data by Reynolds (1997). No clear evidence of warm absorp- 
tion was established (Reynolds, 1997; George et al. 1998), but 
ROSAT data tentatively suggested spectral complexity in the soft 
X-ray (Turner et al. 1993b). Nandra & Pounds (1994) found ev- 
idence of a hard X-ray edge in the Ginga data. The presence 
of the Compton reflection component and a strong Fe Ka line 
was revealed in BeppoSAX data (Bianchi et al. 2001) and later 
confirmed in XMM-Newton observations (Gondoin et al. 2002, 
Longinotti et al. 2007, O'Neill et al. 2007). In 2006 the source 
was observed by the Suzaku satellite (Larsson et al. 2007). The 
spectrum was characterised by a strong soft excess and a reflec- 
tion component, closely resembling the XMM-Newton data; but 
despite the availability of the high-energy data up to 40 keV, it 
has not been possible to model the broadband spectrum unam- 
biguously. In May 2007, Swift caught Mrk 335 in a historical low 
X-ray flux state, with an extremely hard X-ray power law above 
2 keV (Grupe et al. 2007). Following this discovery, the source 
was observed again by XMM-Newton in July 2007 as a target 
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Table 1. Observation Log of Mrk 335 as observed by XMM-Newton 



OB SID 


Date 


RGS exp 


pn exp 


r * soft 


FlllXo.3-2 




Flux 2 _io 




(dd/mm/yyyy) 


(ks) 


(ks) 




(10~ 12 ergs crrr 2 s~') 




(10 -12 ergs cm -2 s~') 


0510010701" 


10/07/2007 


22 


15 


2.84±0.03 


1.95±0.04 


1.02±0.07 


-i ^4+0.23 
-0.10 


0306870101* 


03/01/2006 


130 


92 


2.72±0.01 


32.16±0.04 


2.09+0.01 


17.72±0.07 


0101040101 c 


25/12/2000 


34 


30 


2.80+0.01 


31.74+0.06 


2.16+0.02 


15.00+0.02 



a) Grupe et al. (2008), b) O'Neill et al. (2007), c) Longinotti et al. (2007). The X-ray properties have been measured from the pn spectra. 



of opportunity (ToO) (Grupe et al. 2008). This paper reports on 
the high-resolution spectrum obtained by the Reflection Grating 
Spectrometer (RGS) onboard XMM-Newton. 

2. Observation and data reduction 

Mrk 335 was observed by XMM-Newton on 10 July 2007 
(OBSID 05 1 00 1 070 1 ) for a duration time of 25 ks. Data from the 
three instruments EPIC (0.3-10keV), RGS (5-38 A), and Optical 
Monitor (OM, 6 filters in the 1800-6000 A bandpass) were avail- 
able (Struder et al. 2001; Den Herder et al. 2001; Mason et al. 
2001). The detailed analysis of EPIC and OM data will be pub- 
lished in Grupe et al. (2008). In this paper, only the RGS data 
from the ToO observation are considered. However, Table [Tire- 
ports the broadband properties from all the three XMM-Newton 
available data sets to allow the reader an immediate comparison 
among the flux states of this source. The raw data were processed 
with SAS 7.1.0 with the task rgsproc for the RGS data, which 
produces spectral products for the source and the background 
and response matrices. No background flares due to high-energy 
particles are present in this observation. For the OM data, we 
decided to extract the magnitudes and fluxes from the pipeline 
products (PPS). 

3. RGS spectral analysis 

The spectral analysis was performed using the fitting package 
SPEX[jVer. 2.0), simultaneously fitting both spectra from the 
two RGS cameras. A preliminary analysis was carried out us- 
ing the first-order spectra, re-binned according to the resolution 
of the instrument and with a signal-to-noise ratio of 5 over the 
whole 5-38 A range. This binning was used only to test the 
agreement between the power-law slopes of the RGS and of the 
EPIC data and to visually check for the presence of broad fea- 
tures in the spectrum (Fig. [TJi. Otherwise, since many spectral 
channels fall in the limit of low number of photons, the unbinned 
spectrum was used for the spectral analysis and the C statistic 
was applied (Cash, 1979). The quoted errors correspond to the 
lcr level (i.e. AC=1, for one interesting parameter). Galactic ab- 
sorption of column density N# =3.9x 10 2() crrr 2 is included in all 
the models (Dickey & Lockman, 1990). 

At first glance, the RGS spectrum of Mrk 335 appears to be 
characterised by a weak continuum and a number of emission 
features (Fig. [TJ). The continuum has been fitted with a power 
law with F=2.75±0.17, in agreement with the slope measured by 
the pn instrument (Table [TJ. All the transitions listed in Table [2] 
have been included in the model to fit the narrow emission lines 
in the spectrum. Each line was fitted with a delta line model, 
fixing the centroid energy to the laboratory values. Lines fluxes 
are reported in Table [2] All but two lines, for which only upper 
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Table 2. Fluxes of the soft X-ray lines found in the RGS spec- 
trum. 



Transition A (lab) Flux' f W hm=o AC Flux 2 fwhm(0 vih) 





(A) 


ph m 2 s 1 






ph m 2 s 1 


FeXVII 


17.073 


<0.05 


AG 


= 1 


0.06±0.04 


OVII HeyS 


18.627 


<0.03 


AC= 


=0 


0.06±0.04 


OVIII Lya 


18.969 


0.20+0.04 


AO 


:31 


0.27±0.05 


OVII Hea (r) 


21.600 


0.10±0.06 


AO 


=4 


0.17+0.08 


OVII Hea (i) 


21.790 


0.26±0.08 


AC= 


:21 


0.31±0.09 


OVII Hea (f) 


22.101 


0.17±0.07 


AC= 


:12 


0.16+0.07 


NVII Lya 


24.781 


0.13±0.05 


AC= 


-12 


0.18+0.05 


CVI Ly/3 


28.446 


0.13±0.06 


AC= 


:13 


0.20+0.07 


CVI Lya 


33.736 


0.15+0.06 


AO 


=8 


0.24+0.07 



1) Line width fixed to FWHM=0; 2) Line width fixed to 
FWHM=2200 km s" 1 . The energy centroids were fixed to the 
laboratory value in both cases. The column with the improvement in 
terms of AC refers to the detection significance for the fit with the 
zero- width lines. 



limits are found, yield an improvement in the fit higher or equal 
to AC=4, corresponding to 95.4% for one free parameter. The 
final fit statistic including all transitions in Table[2]is C-statistic= 
6247 for 5264 degrees of freedom (d.o.f.). As a second step, we 
checked for the presence of the radiative recombination continua 
(RRC) from CV, NVII, OVII, and OVIII at ,1=31.622, 18.587, 
16.771, and 14.228 A, which are visible in the binned spectrum. 
They were incorporated in the model, but when this is fitted, only 
upper limits on the emission measures could be found, so they 
are not considered in the following. A careful visual inspection 
of the spectrum did not reveal any absorption feature. 

The detected lines were then checked for line broadening. 
The OVIII Lye line is the most intense isolated transition in the 
spectrum, and it is not affected by bad pixels. A Gaussian pro- 
file was fitted to this line with the wavelength frozen to the ex- 
pected value of 18.96 A and the width free to vary. The measured 
line FWHM is 0.14+0.05 A. We fitted all the spectral lines fix- 
ing their FWHM to the value corresponding to the one of the 
OVIII Lya, obtaining the fluxes reported in the third column of 
Table [2] The spectrum was also checked for line shifts by leav- 
ing the centroids of the strongest transitions free to vary. For 
the OVIII Lya line, a blueshift of A^=-0.014+0.013 A is found, 
with an improvement of AC =7. For the OVII intercombination 
line the 90% limits on the shift with respect to the laboratory 
wavelength are +0.018 and -0.054 A. Taking into account that 
the present signal-to-noise ratio does not allow us to distinguish 
line contamination and that the RGS systematic error of 8 mA is 
included in the error estimate, we conclude that the peak wave- 
lengths are consistent with the laboratory values. 
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Fig. 1. The RGS spectrum of Mrk 335 (corrected for the source redshift). The fitted model, plotted as a red line, consists of a power 
law with F~2.7 and the emission features reported in Table [2] The data have been binned according to the instrumental resolution 
only for plotting purposes (see text). 



3.1. A close look at the H-like and He-like emission lines in 
the spectrum 

The spectral fit in the previous section and the fluxes reported 
in Table [2] reveal that the OVII triplet is characterised by dom- 
inant intercombination line emission and weaker resonant and 
forbidden lines (see Fig.[2]for a zoom on this spectral portion). 
The line ratios in He-like ions provide a powerful diagnostic of 
the physical properties of the emitting gas, mainly by indicating 
whether the X-ray lines are produced by collisional plasma, by 
photoionisation, or in a hybrid environment (Porquet & Dubau 
2000; Porter & Ferland 2007). The lines ratios from the present 
spectra were calculated with the OVIII Lya- and the OVII triplet 

OVIII Lya 

= 1.17 + 0.53, 



lines. The numbers estimated are 



G: 



(/ + 



4.30+2.70, L= 



OVII(f) 
=0.38+0.25, and R= 



:0.65+0.33, 



r 1 1 

where r,/,/ indicate the intensity of the resonant, intercombina- 
tion, and forbidden line components in the triplet. Given the high 
value of the G ratio, albeit measured here with considerable un- 
certainties, the reported line ratios are consistent with an origin 
in photoionised gas and a possible contribution from collisional 
processes (Porquet & Dubau 2000). 

Extensive simulations were run through the photoionisation 
code CLOUDY (Ferland et al. 1998) in order to constrain the 
physical properties of the emitting gas. The aim of these sim- 
ulations was to extract the expected theoretical values for the 
line fluxes and to compare them to the observed fluxes. In addi- 
tion to the line ratios cited above, we also consider the following 
CVILya nric , nnn CVILy/3 



ratios: 

OVIIILya 
NVIILya 



= 0.75 + 0.33, 



OVIIILya 



= 0.65 ± 0.33, 



= 0.65 + 0.28, and 



CVILy/3 



= 0.86 + 0.53. These 



OVIIILya CVILya 
line ratios were calculated assuming the fluxes from the first col- 
umn in Table|2]i.e. assuming lines with a delta profile. The fluxes 



* 1.5x10 - 




Fig. 2. Zoom of the RGS spectrum in the energy range contain- 
ing the emission lines from oxygen. 



measured with non zero-width Gaussian profiles are nonetheless 
consistent with this estimate. 

The theoretical line fluxes were calculated assuming the 
spectral energy distribution (SED) of a standard AGN as cal- 
culated by Korista et al. (1997). The treatment of the UV ra- 
diation field is particularly important when dealing with a high 
ratio between the intercombination and the forbidden lines. In 
fact, an intense UV source can provide enough photons to 
depopulate the forbidden level 2 3 S via photo-excitation and 
then pump the electrons into the intercombination level 2 3 P, 
resulting in a more intense intercombination line (Mewe & 
Schrijver 1978; Kahn et al. 2001). We therefore compared the 
observed SED of Mrk 335 to the one assumed by CLOUDY. 
The UV spectrum as observed by the XMM-Newton Optical 
Monitor is in good agreement with the SED by Korista et al. 
(1997). The spectral slope of the UV-X-ray power-law com- 
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U = 0.8 _ 




log(n„) 

Fig. 3. Contour plots for electron and column density of the line 
emitting gas. The contours were obtained by comparing the ra- 
tios of the emission lines listed in Table[2]to the values predicted 
by CLOUDY simulations (see text for details). The curves define 
the 90 and 99% regions in the parameters space for for all the 
eight line ratios. Only those values of the ionisation parameter 
yielding a valid solution have been plotted. The corresponding 
physical quantities are listed in Table [3] 



Table 3. Summary of the physical properties of the gas clouds 
estimated from the 99% contour plot in Fig. [3] 



logU 


n e range 


rin range 


dist range 


size range 




crrT 3 


-2 

cm 


cm 


cm 


0.4 


10 9.8-1U 


JQ19.9-21.9 


1.95X10 16 - 17 


10 8-12 


0.8 


10 9.5-11 


jq20.3-22 


1.95X10 16 - 17 


10 9-13 



lutions distributed respectively at 90 and 99% considering the 
eight line ratios simultaneously. 

The measured data points are consistent only with CLOUDY 
solutions in the range logU=0.4-0.8, corresponding to the con- 
tinuous and dashed lines in Fig. [3] Since one of the output of 
CLOUDY is the column density of the gas, an estimate of the 
size of the emitting material I was provided using the relation 
l=N H /n e . For each solution, the distance ro of the cloud from 
the ionising source can be extracted from CLOUDY by replacing 
the AGN bolometric luminosity in the definiton of the ionisa- 
tion parameter given above. The value L/, o /=10 447 ergs s -1 was 
adopted after Woo & Urry (2002). For each ionisation parameter 
U, the values of the distance of the emitting gas cloud(s) from 
the nucleus and its size have been estimated. The minimum and 
maximum values correspond to the 99% contour for the range in 
electron and column density of the gas in Fig. [3] The solutions 
admitted by our data are summarised in Tableland they will be 
discussed in the following section. 



F v (2keV) 



monly defined as a ox = -0.3838 log — — — was esti- 

1^(2500 A). 

mated using the pn and OM fluxes at 2 keV and at 2340 A, 
respectively giving F2jtev=0.53xl0~ n ergs crrT 2 s _1 keV -1 and 
F 23 4o=2.5xlO- 14 ergs cirT 2 s _1 A -1 . The resulting UV-X-ray 
power law a ox =-l.32 is fully consistent with the average Seyfert 
value assumed by CLOUDY, i.e. a ox ~\.A. 

The simulations were carried out by varyingthe column den- 
sity of the gas N# and ionisation parameter UQ over a range of 
electronic density n e spanning 10 6 to 10 14 cm 4 . These quan- 
tities were varied in small steps, so that a fine grid of possible 
solutions was obtained. In particular, the ionisation parameter is 
stepped between logU=-2 and log U=2. 

To compare the measured line ratios to the spectral model 
simulated by CLOUDY, each line ratio was searched over the grid 
of predicted values in the following way. For each line ratio i the 
observed value L , with its error err(L o l ) and the corresponding 
theoretical value L t j have been taken into account and used to 

calculate the quantity <x 2 =^— ^ ^r— 

err(L oJ Y 

Then, the quantity cr 2 ) 1 ^ 2 has been minimised, i being 
the summation index over the considered line ratio. The con- 
tour plots in the N/,-n c plane have been drawn by imposing a 
distance equal to 4.61 and 9.21 from the minimum solution, cor- 
responding to 90 and 99% levels of confidence. The inner and 
outer contours in Fig. [3] therefore represent the locii of valid so- 



U is defined as 



<P(H) 



n e c 



where the radiation field is expressed by 



<P(H) 



Anr\ 



Jvl 



dv with vl = lRyd and v2=oo and the relation 



between the normalisation costant fc, and the bolometric luminosity of 
the source is expressed by L = k I nF v dv. 



4. Discussion 

The RGS spectrum of Mrk 335 at low state provides a unique op- 
portunity for getting insights into an X-ray reprocessing region 
that is not easily accessible in Seyfert 1 objects. The soft X-ray 
emission lines could be observed because of the decrease in the 
continuum flux. The nuclear power was attenuated for reasons 
as yet unknown, however the most straightforward explanation 
could be partial obscuration of the central source. At first glance, 
this object seems to be analogous to the other well-known case 
of an extremely variable Seyfert 1, NGC 4051, which showed a 
very rich emission line spectrum when it was observed in a low 
flux state by XMM-Newton (Pounds et al. 2004). According to 
these authors, the analysis of the emission lines in this source, 
especially of the OVII triplet, pointed to interpret the soft X-ray 
spectral features as arising from photoionised low-density gas 
distributed on a large scale, consistent with the AGN narrow line 
region. The variability history of Mrk 335 and the comparison 
with other variable AGN are not the prime focus of the present 
paper. However, the similarity between NGC 4051 and Mrk 335 
suggests to consider first the geometry of the system. 

Mrk 335 is a pure Seyfert 1 galaxy, since the nuclear emis- 
sion does not suffer any obscuration nor does absorption when 
the source is observed at high state, as shown by previous, recent 
X-ray observations (Longinotti et al. 2007, O'Neill et al. 2007). 
Nonetheless, the low-state, soft X-ray spectrum resembles the 
RGS spectra recently observed in obscured AGN by Guainazzi 
& Bianchi (2007). It is then reasonable to postulate that we are 
observing the same line-emitting gas in Mrk 335 as became ob- 
servable due to the drop in the continuum flux. There are some 
objections to this hypothesis. The main point is that, in many ob- 
scured sources, the soft X-ray photoionised matter was demon- 
strated to be coincident with the extended narrow line region, i.e. 
well out of the obscuring torus (Bianchi et al., 2006). The scale 
of this gas is much larger than the distance of the emitter ob- 
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Fig. 4. The cartoon shows qualitatively the location of the 
ionised gas in the nucleus of Mrk 335 on a scale of gravitational 
radii and the corresponding cm scale. 



served in Mrk 335 for all the possible solutions listed in Table[3] 
For this reason, the hypothesis of a common location of the soft 
X-ray photoionised gas for unobscured and obscured objects is 
discarded in the present case. 

The simulations run with CLOUDY cannot uniquely con- 
strain the ionisation parameter or the electron density of the 
gas, and yet they provide a range of physical conditions for the 
gas that are noticeably in agreement. The distance of the pho- 
toionised plasma from the nuclear source ranges approximately 
from 7 to 77 light days, meaning that the gas is always con- 
fined within ~0.06 pc. The column density never gets higher 
than 10 22 cirT 2 (see contours in Fig. [3) and the size of the 
emitting region(s) always appears quite small, the dimension 
of the largest cloud being ~10 n cm. From the black hole mass 
M=14.2+3.7xl0 6 MO (Peterson et al. 2004), the Schwarzschild 
radius R s is estimated to be around 4xl0 12 cm. The distance of 
the Broad Line Region (BLR) from the central source can be 
found by considering the width of the optical lines and assum- 
ing that the BLR gas is virialised. The H/3 FWHM was measured 
by Boroson & Green (1992), and it is equal to 1640 km s _1 , in 
agreement with the reverberation mapping measurements pro- 
vided by (Peterson et al. 2004). The BLR can then be located at 
~7xl0 16 cm, and that this gas in Mrk 335 is constituted of sev- 
eral clouds of emitting material forming a clumpy medium was 
suggested by Arav et al. (1997) on the basis of observational 
constraints. 

It is very likely that the X-ray photoionised gas is inner to 
the BLR clouds and it is also reasonable that the two materials 
form a continuous distribution of clouds. One can imagine that 
the innermost clouds spread out towards the nucleus, reaching 
an ionisation level high enough to emit the observed soft X-ray 
lines. In this case, the spectral line should suffer line broadening 
due to the vicinity of the central black hole, although this effect 
should not be extreme because the Schwarzschild radius is 3 or- 
ders of magnitude smaller. The Gaussian line width measured in 
the OVIII Lya corresponds to a FWHM of 2200+750 km s^ 1 . We 
checked for the width in the OVII intercombination line, and the 
upper limit on the line width is consistent with the measurement 
from the OVIII Lya at <0. 15 A. The distance estimated from the 
virial assumption is then constrained within ~2.3xl0 16 cm and 
1.2xl0 17 cm. In principle, this value can be easily reconciled 
with the distance inferred from CLOUDY; on the other hand, it 
may add an independent observational constraint on the structure 
of the gas. In fact, the underlying assumption in our CLOUDY 
model is that all the X-ray emission lines originate in the same 
cloud, therefore the distances in Table[3]must refer to the bulk of 
the photoionised gas. The distance estimated from the OVIII Lya 
width may indicate that the spectral lines come from a distribu- 
tion of clouds and that OVIII is only concentrated in the inner- 



most ones. To conclude, the X-ray photoionised gas is likely to 
be located within the optical BLR, as depicted qualitatively in 
Fig. [4] Longer exposures of this AGN at a serendipitous low 
state are needed to reach firm conclusions. 

Finally, it has to be remarked that the existence of circum- 
nuclear ionised gas in the innermost region is not a novelty in 
this AGN. Longinotti et al. (2007) reported on the detection of a 
narrow absorption feature around 5.9 keV in the EPIC spectrum 
of the 2000 observation (see Table[TJ. This feature was identified 
as a redshifted Fe XXVI Ka transition and interpreted as the sig- 
nature of gas inflowing towards the nucleus at a velocity of 0.1 1- 
0.15 c. Six years later, XMM-Newton observed the source again; 
the detection of an emission line at ~1 keV in the EPIC spectrum 
was associated to the presence of highly ionised material on a 
distance scale much larger than those discussed herein (O'Neill 
et al. 2007). Interestingly, the most recent X-ray observations of 
Mrk 335 performed during the past year by Suzaku and XMM- 
Newton did not show any emission line at 7 keV (Larsson et 
al. 2007, Grupe et al. 2008), implying that rapid line variability 
must be involved. The potential physical connection among the 
briefly outlined observational results and the soft X-ray lines in 
the present work should be a point of major interest in the next 
studies of this intriguing active galaxy. 
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